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ABSTRACT 

We present an IDL graphical user interface-driven software package designed for the analysis of 
extrasolar planet transit light curves. The Transit Analysis Package (TAP) software uses Markov 



Chain Mon te Carlo (MCMC) techniques to fit light curves using the analytic model of Mandel and 



Agol ( 2002 1. The package incorporates a wavelet based likelihood function developed by Carter and 



Winn (|2009|) which allows the MCMC to assess parameter uncertainties more robustly than classic \ z 



methods by parameterizing uncorrelated "white" and correlated "red" noise. The software is able to 
simultaneously analyze multiple transits observed in different conditions (instrument, filter, weather, 
etc) . The graphical interface allows for the simple execution and interpretation of Bayesian MCMC 
analysis tailored to a user's specific data set and has been thoroughly tested on ground-based and 
Kepler photometry. AutoKep provides a similar GUI for the preparation of Kepler MAST archive 
data for analysis by TAP or any other analysis software. This paper describes the software release 
and provides instructions for its use. 
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1. INITIAL CODE RELEASE 

TAP and autoKep are available for download at http://ifa.hawaii.edu/users/zgazaLk/IfA/TAP.html along 
with descriptions and example videos. A full paper will be made available shortly in which we describe comparison 
tests with published ground-based and Kepler light curves, as well as new data. A preprint is available at the TAP 
website. 

In [J3]we provide a manual for the use of the MCMC fitting feature of TAP, and [j4j provides a manual for autoKep, 
which is designed to streamline the reduction of pipeline-reduced MAST Kepler light curves into a form usable by 
TAP or any transit analysis code. 

2. ATTRIBUTION 

If you use TAP in a publication, please cite this paper and the following, which provide direct contributions to the 
code and/or methods used by TAP: 

1. |Mandel and Agol||2002| for the Mandel & Agol transit light curve model. 

2. |Carter and Winn||2009[ for the wavelet-based red noise treatment 

3. |Eastman et~aT|[20TTl for EXOFAST implementation of the Mandel & Agol model 
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3. TAP MANUAL: 

TAP is designed to input transit light curves and fit them using the Mandel & Agol model. The input light curve 
data must be in ascii files in which the first column contains time in Julian Days (HJD) and the second column contains 
flux (normalized to an "out of transit" value of unity). Any extra columns are ignored by TAP. 



1. Starting TAP: 

a) Start an IDL session: 

b) Run 'tap' at the idl command line: This creates an instance of the TAP GUI. No further interaction 

with the IDL command line is necessary. 

2. Load one or multiple transits: 

a) Click on "Transit File": under the "Load Transit" menu 

b) Select your file, click "ok" in the pop-up window: Files must be in ascii format. TAP uses the first 

two columns. First column must be a timestamp in days. Second column must be flux and should be 
normalized to an out of transit value of 1. Additional columns are not used. 

c) Set Resampling mode: For long integration observations (such as Kepler's 29.4244 minutes per datapoint), 

the Mandel and Agol light curve must be calculated on a finer resampled time scale and then rebinned to 
the observed cadence. If this is the case, select the 'Rebin to "Input Integration"' option. In this case, 
the Input Integration box will become active. Set your integration time and the number of samples (per 



datapoint) for resampling. In general, N=4 samples is sufficient for Kepler light curves. See Kipping (20101 
for a more complete description. 

d) Click "Load Transit": The plot window will display the transit light curve and residuals in the bottom 
plot. The initial transit guess will appear as a blue curve, and the residuals will be plotted as distance from 
that blue curve. The default strength for "correlated red noise" is (i.e. no red noise). If you choose to fit 
for red noise as well, a red curve will also be overplotted, showing instead the transit model with the red 
noise signal added to it. 

NOTE: To load multiple transits, either repeat the above steps with each ascii file or prepare an ascii file with 
all transits separated by a line of -l's matching the number of columns in your file. TAP will automatically 
load each transit separately. 

Restore an existing setup: If you have started an MCMC calculation and want to restore that exact setup, 
click on the "Setup File" button under the "Load Transit" menu. Select a file, which TAP automatically 
creates, called "TAP_setup.idlsav" which is automatically saved into the directory of an MCMC run. The 
setup will be restored, and changes can then be made before beginning a new MCMC calculation. 

3. Adjust parameters for multiple transits: 

a) Click on the "Manage Transits" menu button: 

b) General items: Sliders which control which transit is "active" and options for plot scaling. The "active 

transit" is the one which is affected by manual changes to the fit parameters. It is labeled with blue arrows 
in the plot and its current parameters arc plotted in the Parameter box to the left of the plot. 

c) Setup Inter- Transit Links: This button launches a new window which contains an entry for each parameter 

of each transit. The values in the boxes represent the current "set" the parameter belongs to. Parameters 
in any row with the same set are locked together, such that those transits have the same parameter values 
and must evolve together in the MCMC calculations. 



4. Fit the transit light curve: 



a) Click on the "Fit" menu button: 

b) "Manual Parameter Adjustment": this button opens a new window with 13 sliders. These represent the 

possible free parameters for MCMC fitting. The values can be adjusted by dragging the sliders or by typing 
in values and pressing "return". The light curve plotted in the main window will update automatically. 

"System Parameters" The parameters of the transiting system from the Mandel & Agol (2002) analytic 
light curve. 

"Quadratic Limb Darkening" Linear and quadratic Limb Darkening Coefficients. 

"Data Specific Parameters" A linear correction for airmass or normalization. Slope and Y-intercept of 
the "out of transit" level to account for any trend in the data. Also Noise parameters: Strength of "White" 
uncorrelated gaussian noise and "Red" (1/f) correlated noise. 



c) "Parameter Limits and Locks" : this button opens a new window with one row for each of the 13 possible 

free parameters. 

"Min/Max": The minimum and maximum of the sliders for "Manual Parameter Adjustment", and, if 

"Apply Limits to Fitting" is selected, the MCMC is not allowed to wander beyond these limits. 

"Lock": Select this to lock a parameter during fitting. This is used for parameters that are unable to be 

determined from the transit light curve, such as period in the case of a single light curve. 

"MCMC Accept Rate": This is the percent of jumps for which the MCMC algorithm will optimize the 

jump betas. The typical value for efficient MCMC convergence is 44default value. 

"Apply Limits to Fitting": Will force the MCMC execution to explore a parameter space limited by 
the "Max/Min" values. 

d) "MCMC Parameters": Set the number of chains and links per chain for the MCMC execution. 

e) "Gaussian Priors" : Penalty terms allowing a parameter to vary in the MCMC analysis but under a gaussian 

penalty term. The term works as a modification to the Likelihood of any particular fit by adding a term 
of the form exp[— (value — Xi) 2 /sigma 2 ], where value and sigma correspond to the set priors and Xi is the 
current MCMC value of that parameter. This is useful when a parameter is theoretically known but should 
be allowed to wander, most commonly the limb darkening parameters. The code will then find an ideal 
solution between the "Value" and mathematically most likely solution. Set a value, sigma, and click the 
"Enable" box to activate this feature. 

f) "Execute Chain": This will open a new window and freeze the base window. The new window will display 

the progress of the MCMC chain analysis and the plot window on the frozen TAP base will cycle through 
the evolving parameter distributions. The GUI will create a new directory to contain the MCMC run in 
the directory from which TAP was started. 

Analyze MCMC Chain(s) 

a) Click on the "MCMC Inference" menu button: 

b) Select options for output: A .tex document is automatically created containing details of the MCMC run 

and tables of results. Beware of creating ascii files of the MCMC chains-they can become quite large. Each 
output ascii file contains in its header a code snippet to load into IDL using 'readcoP 

c) Click the "MCMC Save File" Button: A dialog will appear to allow you to select a 'TAP .setup. idlsav' 

file. Select this file in any MCMC directory which has completed at least 1 MCMC chain and press "ok". 

d) Click the "Load" button: TAP does the rest, automatically loading, combining, and analyzing the com- 

pleted chains. You will find the output files in the root MCMC directory. You can perform the MCMC 
inference as often as you would like, and as new chains complete the software will automatically utilize the 
full available results. 



■5 



\ TAP MCMC v2. 02 



ansit Analysis Package v2.02 

Quit 




(10;51;ie) Running Chain 03 of 10 

(10:51:16) Stabilizing Betas 
(10:51:32) Coarse Stabilization Done 
(10:52:35) Fine Stabilization Done 

1(10:52:35) Betas have stabilized 



Sigma White: 



0.00100000 



Load Transit | 


Manage Transits | 


| Fit 


MCMC Inference | 



Manual Parameter Adjustments and Setup 


Markov Chain Monte Carlo 


Manual Parameter Adjustment | 




MCMC Parameters 








Adjust Limits and Locks | 




Gaussian Priors 






Execute Chain 




Fig. 1. — The Transit Analysis Package (TAP) Graphical User Interface while executing a Markov Chain Monte Carlo analysis on three 
light curves of WASP-lOb. The main background widget freezes during MCMC calculation as a new smaller widget is spawned in the 
foreground. The MCMC widget provides updates on the status of the analysis process while plotting every tenth x n model state with a 
blue analytic model and red curve representing the model and red noise. The top and bottom progress bar represent the current chain 
and overall completeness, respectively. During MCMC execution, the window in the main widget cycles through the currently explored 
probability density distributions for different parameters so that a user may monitor the progress of the analysis. 
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4. AUTOKEP MANUAL: 

autoKep is designed to load in *llc.fits (pipeline reduced light curve) files from the Kepler MAST archive. The 
transits are output as ascii files and (optionally) loaded into TAP. 

1. Starting autoKep: 

a) Start an IDL session: 

b) Run 'autokep' at the idl command line: This creates an instance of the autoKep GUI. No further in- 

teraction with the IDL command line is necessary. 

2. Load one or multiple Kepler MAST .fits files: 

a) Click on "Fits File" under the "Load Fits" menu: in the pop up dialog, select a Kepler MAST .fits 

file, then press "ok" 

b) Click "Load Fits File" or "Load all .fits": The former option will load the currently selected fits, the 

latter loads all .fits files in the directory from which autoKep was started. 

3. Detect Transits 

a) Click on "Detect Transits" menu item: The GUI workspace will change. 

b) Find the first transit: Either click on "Auto Detect" and give the software a shot at automatically finding 

the transit or click "Manual Detect", then click to the left and right of a transit in the top (white back- 
grounded) plot. In either case, if autoKep finds a plot a popup window will ask you to verify the detection. 
Do so. 

c) Find the second transit: Repeat the procedure by clicking a button in the "Second Transit" submenu. 

You should select an epoch immediately neighboring the initial transit detection, as autoKep will then use 
these two transits to estimate orbital period and "map" the remaining transits in yellow. 

d) Map the transits: If the yellow regions do not cover the transits in the top plots, click on the "Adj. Yellow 

Maps" button. A widget will open which allows you to modify period and duration (and thus the location 
and witdth, respectively, of the predicted transit epochs) anchored to the second detected transit. Once 
the yellow maps look OK, click on the "Map Transits" button. This button searches each yellow "epoch" 
for a transit signature and stores each as a transit light curve. 

4. Adjust the Transits 

a) Click on "Adjust Transits" menu item: The GUI workspace will change. 

b) Adjust the Duration of the Transits: The left two horizontal sliders adjust the duration of the transit. 

The top slider is the transit event itself (red points) . These points are ignored when out of transit variations 
are fit, so it is important to make sure the entire transit is covered. The bottom slider marks the extent of 
the out of transit light curve (blue points) which will be fit with a polynomial trend. Enough points are 
needed to properly fit the curve. All transits are fit with the same duration and OOT windows. 

c) Active transit: the vertical slider sets the active transit, which is plotted in the lower left plot and in color 

in the lower right plot. This slider is useful for assuring that the out of transit trend fits each transit well. 

d) Fit OOT Trends: the rightmost pair of horizontal sliders adjust the polynomial order of the OOT trend. 

The top slider adjusts only the trend of the active transit, while the bottom slider adjusts all transits at 
once (to the same value). 

5. Export the Transits 

a) Click on "Export Transits" menu item: The GUI workspace will change. 

b) Set the "Write Region": The extent (in days) of the exported transit. To reduce MCMC calculation time 

it may be pertinent to cut down the output, being certain to leave sufficient out of transit points for the 
MCMC model fitting. 

c) Export Options: A user may choose whether or not to apply the OOT corrections set in the "Adjust 

Transits" menu, and may also choose to automatically export the transit to TAP for analysis (the file is 
also written as a .ascii file). 

d) Root Filename: Type in the base filename, to which .ascii will be appended, for the output transit. A single 

file is written, with 3 columns: BJD, normalized flux, and photometric error. Each transit is separated by 
a line of -l's, as is the standard form for loading multiple transits into TAP from a single file. 

e) Export Transits: Click the lower Export Transits button. If you have chosen the "Export to TAP" option, 

autoKep will automatically exit and TAP will start. Otherwise, the GUI remains active for additional 
changes. 



autoKep Kepler LC Prep vO.5 
Quit 




A 
J 



(13:32:13) Transit located* Search for 2nd transit* 
(13:32:14) Click before and after a transit ONE epoch away 
1(13:32:16) ...possible transit signature detected 



1,0013 
0.9966 
0.9919 



1 



2 



3 



4 



5 



I 



\ 
t 7 



9 



0.54 



9.43 



18.32 27.21 
Days since 54953.000 



36.09 



44. 




R 1.015 
e 
1 

p 1.010 

1 

u 

X 1.005 



* »• « 



n 0.995 
t 



28.78 29.11 
Days since 54953.000 



29.43 



-10 



-5 5 10 

Hours from Mid Transit 



1) Load Fits | 2) Detect Transits | | 3) Adjust Transits 4) Export Transits | 



0.29727 



I ~ 



Duration 



1.0135 



00T 



Transit 



Fit DOT Trends 
actiue order: 
5 



all order: 
5 



-A 



Fig. 2. — The autoKep Graphical User Interface. 
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